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SUMMARY 

A  series  of  significant  multiwavelength  laser  propagation  experi¬ 
ments  have  been  conducted  with  point-source  transmitters  and  independent 
determinations  of  turbulence  strengths.  Under  non-saturated  scintillation 
conditions,  we  find  that  the  theoretical  wavelength  prediction  of  log  ampli¬ 
tude  variance  is  valid.  Also,  we  find  that  the  near-  and  far-infrared  wave¬ 
lengths  do  not  saturate  at  lower  variances  than  visible  wavelengths,  and 
large  infrared  variances  have  been  measured.  We  find  that  different  tech¬ 
niques  for  measuring  the  strength  of  turbulence  yield  different  (but  related) 
results . 

We  have  found  that  the  multiwavelength  covariance  characteristics 
are  usually  independent  of  conditions  and  are  near  theoretical  values,  although 
certain  anomalous  results  suggest  limitations  in  the  extent  of  the  inertial 
subrange  turbulence  model.  We  have  found  that  large  receiver  apertures 
do  not  result  in  the  degree  of  reduction  in  total  signal  fluctuations  that 
would  be  expected  from  the  theory.  Finally,  we  have  found  that  simultaneous 
records  of  changing  variances  and  turbulence  strengths  yield  interesting 
characteristics  which  may  be  related  to  the  stationarity  of  the  medium. 
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I. 


INTRODUCTION  AND'  SUMMARY 


W 


During  the  period  since  the  preceding  Quarterly  Progress  Report, 
we  have  conducted  a  series  of  significant  experiments  in  multiwavelength 
laser  propagation  with  a  virtual  point-sour ce  (sp?nerical  wave)  transmitter 
configuration.  As  discussed  in  the  preceding  report,  the  customary  use  of 
a  finite,  divergent  transmitter  beam  does  not  result  in  data  which  may  be 
simply  compared  to  the  predictions  of  spherical  wave  theory.  Also,  an 
important  feature  of  the  present  data  is  that  they  are  backed  up  by  several 
independent  determinations  of  the  strength  of  the  turbulence.  Hence,  although 
much  more  data  must  be  gathered  before  firm,  general  conclusions  can  be 
drawn,  the  results  of  our  present  experiments  suggest  new  and  important 
hypotheses . 

In  particular,  during  unsaturated  conditions  we  find  that  the  theoretical 
k7/6  dependence  of  log  amplitude  variance  is  reasonably  well  verified  within 
experimental  spread.  Also,  we  find  that  the  near-  and  far -infrared  wavelengths 
do  not  saturate  at  lower  variances  than  visible  wavelengths,  and  we  have  seen 
variances  as  high  as  0,63  and  0.  15  at  1.15  and  10.6  micron  wavelengths  , 
respectively.  We  find  that  different  techniques  for  measuring  the  strength  of 
turbulence  yield  different  (but  related)  results. 

We  have  found  that  the  multiwavelength  covariance  or  lateral  correla¬ 
tions  of  log  amplitude  are  usually  independent  of  conditions  and  are  near 
theoretical  values,  although  certain  anomalous  results  have  been  noted  as 
discussed  later.  These  anomalies  may  relate  to  a  relatively  large  inner 
scale  (limited  inertial  subrange).  Also,  we  have  confirmed  other  reports 
that  large  receiver  apertures  do  not  result  in  the  expected  degree  of  reduc¬ 
tion  in  total  signal  variance,  and  we  have  therefore  shown  that  this  suprising 
phenomenon  is  not  fundamentally  due  to  the  use  of  a  finite  transmitter  aperture. 
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Finally,  we  find  that  simultaneous  records  of  variance  and  turbulence 
strengths  vs,  time  yield  interesting  characteristics,  which  may  be  related 
to  considerations  of  stationarity  of  the  medium. 

In  the  following  sections,  we  discuss  the  above  aspects  in  detail. 

H .  EXPERIMENTAL  RESULTS 

A.  GENERAL  MEASUREMENTS  PROGRAM 

A  Standard  Operating  Procedure  was  established  in  which  we 
measured  the  following  quantities  during  experimental  "run":* 

1.  Strength  of  turbulence  (O^2)  through  determination 

of  the  mean-square  temperature  differential  between 
two  thermal  microprobes  (<A  Tz>) 

2.  C^2  from  the  slope  of  the  probability  curve  for  AT 
on  gaussian  probability  axes 

3.  C^2  ("optical")  inferred  from  the  variance  (probability 
slope)  of  the  log  amplitude  scintillations  of  a  6328A 
portable  laser  located  at  500  feet  from  the  receiver 

4.  Simultaneous  log  amplitude  probability  distributions 
and  variances  [C^(0)]  over  a  one  mile  path,  for 
4880A,  1.15  microns,  and  10.6  microns 

5.  Simultaneous  log  amplitude  covariances  [C^(r)]  for 
the  same  three  wavelengths 
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6.  Simultaneous  speccra  of  amplitude  scintillations 
for  the  same  three  wavelengths 

Items  #1-3  above  are  performed  simultaneously,  followed 
by  #'s  4  and  5  in  quick  sequence,  followed  typically  by  a  repeat  of  1-3.  The 
entire  data  run  takes  less  than  30  minutes.  Wind  speed  and  direction  (at 
both  ends  of  the  path),  temperature,  and  barometric  pressure  are  recorded 
continuously . 


The  inverse  Fresnel  number  for  the  4880A  transmitter  is  on 
the  order  of  104  and  for  10.6  microns  is  102.  At  1.15  microns,  the  transmit¬ 
ter  is  larger  but  contains  several  transverse  modes  and  hence  may  be  expected 
to  yield  results  which  are  in  good  agreement  with  those  from  a  point  source; 
this  was  borne  out  by  the  experiments.  The  receiver  aperture  at  all  wave¬ 
lengths  is  3mm. 

Preliminary  experiments  were  also  performed  for  ancillary 
quantities,  as  follows: 

1.  Log  amplitude  probability1  distribution  and  variance 
at  4880A  vs.  receiver  aperture  up  to  121/,;  inches 
diameter  . 

2.  High-frequency  cut-off  of  the  spectrum  of  tempera¬ 
ture  fluctuations,  as  an  indication  of  the  size  of  the 
inner  scale 


❖  This  determination  is  not  presently  feasible  for  very  light  transverse  wind 
conditions,  since  our  spectrum  analyzer  is  inaccurate  at  or  below  10  Hz. 
However,  we  are  obtaining  an  analyzer  of  greater  low-frequency  resolution. 
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3. 


Simultaneous  recordings  of  the  log  amplitude 

variances  and  C  2  vs.  time, 
n 

As  appropriate,  these  measurements  will  be  incorporated  in  I 

the  Standard  Run  described  above.  Also,  at  a  later  date,  transmitter  apertures 
will  be  systematically  increased.  However,  at  the  present  time  we  are  con¬ 
centrating  on  obtaining  a  sufficient  quantity  of  data  to  support  preliminary 
conclusions  . 


B.  NATURE  OF  RESULTS 

A  typical  oscillogram  of  three -wavelength  amplitudes  is  shown 

in  Figure  1.  Unlike  earlier  results  for  finite  apertures,  the  point-transmitter 

scintillations  shown  here  are  only  weakly  correlated  between  wavelengths. 

Except  for  highly  non -stationary  conditions,  the  log  amplitudes  were  always 

log  normal(for  both  saturated  and  nonsaturated  conditions )  within  the  accurate 

probability  resolution  range  of  our  technique  (Fig.  2).  We  can  expand  the 
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probability  tails  at  a  later  date,  should  this  prove  desirable. 

Covariance  results  were  relatively  consistent  and  independent 
of  conditions,  including  the  saturation  of  scintillations.  We  define  the  trans¬ 
verse  leg  amplitude  correlation  length  r  as  that  (horizontal)  separation  of 
detectors  which  results  in  the  following  condition: 

c,<ra>=  ^=,(0)  (1) 

The  experimental  value  of  r  was  generally  somewhat  smaller  than  that  of 
the  theory  of  Ref.  3,  especially  at  10.6  microns.  A  typical  covariance  result 
is  shown  in  Figure  3.  Major  anomalies  in  covariance  were  occasionally 
noted,  and  will  be  discussed  below. 
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The  spectra  of  scintillations  are  exemplified  by  the  smoothed 
curves  of  Figure  4.  Although  the  curves  are  not  readily  normalized  instru- 
mentally,  due  to  the  extrapolation  necessary  to  obtain  the  O-Hz  level,  the 
"bandwidth"  corresponding  to  (1/e)  of  this  level  is  readily  obtained  and  in  all 
cases  agreed  rather  well  with  the  transverse  wind  component  and  measured 
value  of  r  .  This  substantiates  the  "frozen-in"  structure  of  the  scintilla- 

cl 

tion  pattern. 


The  three  methods  of  obtaining  C  2  proved  highly  interesting. 

The  approach  was  suggested  by  the  work  of  Ref.  4.  We  confirm  the  pre¬ 
liminary  findings  of  that  report,  m  which  it  was  stated  that  the  temperature 

$ 

fluctuations  (AT)  depart  seriously  from  a  gaussian  distribution.  Representa¬ 
tive  results  for  three  different  runs  are  shown  in  Figure  5,  showing  (l)  a  good 
gaussian  characteristic;  (2)  a  pronounced  "tail"  all  but  obliterating  a  gaussian 
characteristic;  and  (3)  a  curve  best  represented  as  two  discrete  gaussians 
due  to  a  sudden  change  in  turbulence  conditions  (e.g.  ,  when  the  sun  emerges 
from  behind  a  cloud).  In  all  cases,  the  optical/IR  scintillations  were  log 
normal . 


The  typical  characteristic  evidenced  a  linear  (i.e.  gaussian) 

region,  with  definite  tails,  and  consequently  the  value  of  C ^  obtained  directly 

from  (AT)2  was  always  larger  than  that  from  the  probability  slope.  In  addition, 

the  C  2  inferred  from  the  short-path  optical  scintillations  always  fell  in  between 
n 

the  other  two  values.  The  amount  of  spread  in  the  three  values  depended  on 
the  degree  of  nongaussian  behavior,  and  the  one -mile  scintillations  (e.g.  at 
10.6  microns  ,  which  was  never  saturated)  were  fairly  well  related  to  the  latter 
two  determinations  of  C^2.  Clearly,  the  interaction  of  the  optical  filter 


^However,  we  did  not  see  clear  "spiking"  in  the  temporal  AT  behavior. 
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function  and  the  turbulence  spectrum  should  be  thoroughly  investigated.  This 
will  be  discussed  more  fully  below. 

C.  REVIEW  OF  RESULTS 

Table  I  contains  a  review  of  results  from  selected  runs  which 
illustrate  various  aspects  of  the  phenomena. 

It  should  be  pointed  out  initially  that  the  accuracy  inherent  in 
the  instrumentation  and  approach  is  very  good,  since  very  large  dynamic 
ranges,  and  probability  slope  methods,  are  used.  However,  it  is  character¬ 
istic  of  atmospheric  phenomena  that  atationarity  cf  the  medium  is  a  useful 
approximation  at  best,  and  hence,  the  spread  in  results  will  always  be  sub¬ 
stantial.  This  will  be  directly  illustrated  in  a  later  section. 

In  Table  I,  the  theoretical  values  of  C^(0)  are  obtained  from 
the  standard,  Rytov-based  theory  for  a  point  source  (Ref.  3}.  Hence,  these 
values  should  not  be  expected  to  predict  saturation,  i.e.,  the  limitation  of 
experimental  C^(0)  to  values  on  the  order  of  0.6. 

Finally,  we  point  out  that  although  our  range  is  nominally  flat 
and  uniform,  it  is  not  as  ideal  in  this  regard  as  e.g.  the  Emerson  Dry  Lake 
site  (Ref.  4).  In  the  case  of  a  Westerly  wind  (Runs  #5  and  6),  there  are  up¬ 
wind  trees  over  th*>  portion  of  the  path  near  the  receiver.  These  character¬ 
istics,  which  are  logistically  necessary,  do  not  comprise  a  serious  limitation 
to  the  goals  of  this  program,  since  the  distribution  of  any  nonuniformity  in 
over  the  path  weights  each  wavelength  identically.  Hence,  although  this 
may  contribute  seme  to  the  apread  in  the  relationship  between  independently 

measured  C  2  values  (at  the  receiver)  and  nonsaturated  scintillation  variances, 
n 

the  validity  of  the  theoretical  k7/6  dependence  of  variance,  and  the  wavelength 
dependencies  of  saturation,  are  investigated  without  difficulty. 
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We  now  point  out  the  salient  implications  of  each  run. 

Run  #1  may  be  characterized  as  "typical"  of  afternoon  data  on 
a  clear  winter  day  with  some  wind.  Scintillations  at  4880A  are  essentially 
saturated,  and  at  other  wavelengths  C^2  from  (AT)2  predicts  variances  of 
2-3  times  more  than  the  observed  values,  due  to  the  nongaussian  behavior  of 
AT.  The  measured  values  of  ra  agree  with  theory  for  4880A  and  are  some¬ 
what  low  at  10.6  p.  The  scintillation  bandwidth  is  reasonably  consistent  with 
r&  and  the  transverse  wind  velocity. 

In  Run  # 2,  there  was  very  little  wind  and  hence  there  may  have 

4 

been  insufficient  energy  for  broad  applicability  of  the  inertial  subrange.  The 
most  noticeable  anomaly  is  the  unusually  high  value  of  variance  at  4880A, 
a  result  which  was  obtained  with  repeated  trials.  The  variances  at  the  infrared 
wavelengths  are  accurately  predicted  from  C^(C)  or  the  "optical  C^2."  This 
value  agrees  within  50%  of  that  determined  from  the  AT  probability  slope, 
while  that  from  (AT)2  is  3.6  times  higher.  Values  of  r  a>’e  typical.  A  further 

ct 

evidence  of  a  large  inner  scale  for  this  run  will  be  discussed  later. 

In  Run  #3,  a  higher  wind  was  sustained,  and  scintillations  were 
not  saturated.  The  ratio  of  variances  is  in  very  good  agreement  with  the 
k7/6  prediction.  However,  there  are  considerable  tails  on  the  temperature 
probability  distribution,  and  a  large  value  for  C^2  (A).  The  variances  agree 
most  closely  with  the  predictions  from  Cn2(C).  Values  of  r&  are  typical,  and 
scintillation  bandwidths  and  transverse  wind  are  consistent  with  these  values. 

Run  #4  is  typical  of  data  being  obtained  as  the  spring  weather 
becomes  warmer,  and  4880A  scintillations  are  saturated.  Very  severe  tails 
on  the  AT  distribution  were  evidenced,  and  there  was  a  large  spread  between 
the  three  measurements  of  Cn2,  with  the  optical  value  fairly  well  predicting  the 
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observed  variances.  The  low  wind  and  high  value  of  r^  at  4880A  may  indicate 
a  large  Inner  scale. 

In  run  #5,  taken  in  the  evening,  the  wind  wai  alight  but  the 
tail  in  AT  and  spread  in  C^2  measurements  were  relatively  small.  4880A  was 
saturated,  and  the  other  variances  were  best  predicted  by  the  optical  C^1,  or 
even  tlu-t  from  (AT)2. 

Finally,  run  #6  was  taken  during  the  early  morning  under 
relatively  light  turbulence  conditions.  The  turbulence  was  considerably  non¬ 
stationary  ,  and  values  of  r  wero  quite  anomalous,  indicating  a  breakdown  in 

ft 

the  subrange.  The  variances  (a)  were  taken  from  probability  curves,  which 
were  somewhat  erratic  due  to  nonstationarlty;  the  values  (b)  were  taken 
directly  from  dynamic  recordings  of  C|(0)  vs.  t  at  coincident  relative  maxima, 
as  discussed  later.  The  latter  values  are  most  closely  related  to  C  *(C). 

The  results  of  Runs  2,4,  and  5  indicate  that  saturation  does  not 
occur  at  lower  variances  for  longer  wavelengths .  This  is  further  supported  by 
a  measured  variance  at  10. 6^  of  0.15,  as  discussed  below.  This  conclusion 
contradicts  highly  conjectural  statements*  made  previously,**10  and  would 
resolve  a  serious  dimensional  difficulty.10’11  The  result,  if  continually 
supported  by  further  experiments ,  ha*  important  implications  for  systems 
operating  at  10.6  microns,  since  the  dynamic  range  of  fading  is  given  approxi¬ 
mately  by 

D.R.  -  100  Ct(0)  dB  (2) 


*In  Ref.  10,  these  statements  were  based  on  sparse  experimental  data  with 
no  C  *  backup.  In  Ref.  4,  the  data  on  IR  scintillations  were  again  sparse, 
and  &e  variances  predicted  from  Cn*  were  not  large  enough  to  test  the  presence 
of  IR  saturation. 
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For  instance,  under  the  conditions  of  runs  #2  and  4,  and  the  high  10,6  micron 
'variance  reported  beiow,  saturation  at  the  t0j  =  0.6]  level  is  distinctly 
possible  for  10.6  micron  propagation  over  a  low,  10  mile  path. 


D. 


SIZE  OF  THE  INNER  SCALE 


. .  712  4,9 

:  It  has  been  recognized  by  theoretical  '  and  experimental 

investigators  that  the  extent  of  the  inertial  subrange  should  be  known  for  each 

experimental  condition.  That  is  ,  although  the  Kolmogorov  turbulence  spectrum 

lead#  to  results  of  obvious  utility  (and  there  exists  no  useful  alternative  model) , 

various  effects  under  certain  diurnal  and  wind  conditions  suggest  limitations 

QH  the  subrange,  which  in  turn  lead  to  anomalies  in  measured  variances  ,  co- 

variances,  etc.  In  particular,  it  would  be  highly  desirable  to  know  the  non- 
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zero  value  of  the  inner  scale,  e.g.  for  use  with  the  optical  filter  function. 

3 

It  may  be  pointed  out  that  the  theoretical  covariance  result  is  valid  only  when 
the  receiver  is  in  the  far  field  of  the  inner  scale. 


Although  the  scintillation  phenomena  may  be  used  to  'probe1 
7 

the  subrange  ,  for  example  using  a  variable  transmitter  aperture,  a  more 

direct  approach  for  our  purposes  would  be  to  determine  the  inner  scale  from 

12 

the  spectrum  of  temperature  fluctuations.  We  are  currently  working  on 
this  problem  experimentally,  using  our  thermal  microprobe,  and  if  successful, 
we  will  inforporate  this  determination  into  the  Standard  Operating  Procedure. 


Am  indication  that  the  method  may  be  successful  may  be  seen 
from  Figure  6.  Although  the  theoretical  rms  spectrum  falls  off  as 
(frequency)  5/&,  a  distinct  flattering  -  out  is  observed  well  above  the  instru¬ 
mental  noise.  This  character istic ,  which  i3  for  Run  #2  above,  implies 


♦An  attempt  to  directly  measure  the  breakdown  of  the  "2/3  law"  dependence 
of  the  mean  square  differential  temperature  between  two  points  is  hampered 
bv  the  impra  cti  cal  ity  of  bringing  probes  suitably  close  together,  especially 
without  disturbing  the  temperature  field. 


* 


\ 


A 
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an  inner  scale  of  approximately  2  cm  when  the  average  wind  velocity  is  ta.'en 
into  account. 


We  will  improve  the  method  by  (1)  utilising  a  spectrum  analyser 
of  greater  capability  at  very  low  frequencies;  and  (2)  utilising  discrete 
spectral  points  over  a  logarithmic  amplitude  and  frequency  scale.  It  is 
necessary  that  the  full  dynamic  range  of  the  spectrum  analyser  be  utilised, 
so  that  AT  must  be  taken  directly  from  the  microprobe  (Fig.  7)  rather  than 
from  an  FM  magnetic  tepe  recording. 

E.  APERTURE  AVERAGING 

13  4 

It  has  been  reported  by  Fried  and  co-workers  ’  that  un¬ 
expectedly  large  fluctuations  remain  when  the  scintillating  beam  is  received 
by  a  large  aperture.  It  was  speculated  that  the  effect  relates  to  finite  transmit¬ 
ter  aperture  sir.es. 

We  have  conducted  preliminary  receiver  aperture  averaging 
experiments  in  order  to  determine  the  presence  and  extent  of  this  effect  for 
our  point  -source  transmitter  at  4880A.  (The  work  will  be  extended  to  other 
wavelengths  at  a  later  date.)  We  find  that  the  aperture-averaged  reduction 
in  total  signal  variance  is  somewhat  better  than  that  reported  in  Ref.  4,  but 
nevertheless  that  a  substantial  fluctuation  remains.  Furthermore,  the  effect 
is  observed  when  simultaneous  measurements  of  covariance  indicate  a  aero 
corrclat.oo  for  r  >  ra ,  so  that  e.g.  beam  wander  cannot  explain  the  phenomenon. 
Thus  the  effect  retrains  a  mystery. 

The  aperture-averaging  results  are  shown  for  two  runs  in 
Figure  8.  The  results  were  independent  of  the  detector  configuration  (focused 
photodiode  vs.  unfocused  PMT).  The  data  for  typical  covariance  conditions 
and  a  saturated  variance  show  a  residual,  for  the  total  aperture  averaged 


variance, on  the  order  of  10%  of  C^{0).  The  log  amplitude  probability  distri- 
iutiniLis  gau&sian  (Fig,  9).  The  other  inonsaturated)  data  are  for  the  anoma- 

•.  - -  loualy  large  covariance  conditions  of  Run  ^6,  and,  although  larger  intermediate 

apertures  are  necessary  to  achieve  averagings  the  percentage  reduction  at 
”  full  aperture  is  better  than  that  for  small  covariances .  This  is  indeed  sur - 
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fey  the  incre&&e  in  variance  at  2.5  cm  over  that  for  a  point  detector .  However, 
'the  general  results  of  this  experiment  were  repeatable. 


In  order  to  compare  the  aperture -averaging  results  with  a 
theoretical  prediction,  we  may  utilize  the  results  of  Ref.  14.  Although  this 
analysis  wa 8  done  for  a  plane  wave  source  under  Rytov  (nonsaturated)  condi¬ 
tions,  the  covariance  characteristics  for  a  point  source  and/or  for  saturated 
conditions  are  not  materially  different.  To  summarize  the  results  of  this 
analysis,  we  define  a  transverse  scintillation  correlation  length  r  ’  =  (4L/k)/£ 
-,2.24  cm  for  our  path  at  48S0A.  Then,  although  Ref.  14  deals  directly  with 
Cj(6)  rather  than  C^(6),  it  may  be  deduced  from  the  results  presented  there 
that 

_C  -  (D/r  ')*  C.l 0)  (3) 

L,  3.  X 

large 
aperture 
diameter  D 


This  result  is  as  intuitively  expected  from  the  independent  addition  of  approxi¬ 
mately  (D/r  ')*  independent,  log  normal  quantities.  For  a  12  inch  or  30  cm 

a 

aperture,  we  would  then  expect  a  variance  of  0.005  times  that  for  a  point 
aperture.  This  is  an  order  of  magnitude  less  than  the  lowest  value  actually 
measured. 
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F.  DYNAMIC  RECORDING  OF  C  2  AND  VARIANCES 

n 

As  discussed  in  the  next  section,  the  turbulent  medium  may 
be  highly  nonstationary.  The  reason  may  be  obvious,  as  in  the  case  of  small 
clouds  periodically  obscuring  the  sun;  or  more  subtle,  as  in  the  case  of 
Run  #6. 


The  three -wavelength  variances  vs.  time  for  Run  #6  are  shown 
in  Figure  10.  The  variances  are  well  correlated,  but  it  is  obvious  that  their 
ratio  is  nonconstant,  even  though  conditions  were  nonsaturated.  This  un¬ 
doubtedly  relates  to  the  changing  extent  of  the  inertial  subrange.  For  multi¬ 
wavelength  comparisons  under  such  unstable  conditions,  such  a  recording 
permits  the  use  of  truly  simultaneous  values,  e.g.  at  a  relative  maximum  or 
minimum. 


A  recording  of  thermal  [  (AT)*  J  Cft*,  optical  C^*  from  a  short 

path,  and  variances  at  4880A  and  10.6  microns  is  shown  in  Figure  11.  For 

this  recording,  the  wind  was  blowing  parallel  to  the  path,  from  the  receiver 

(where  the  C  *  instrumentation  is  located)  to  the  transmitter,  and  the  sky 
n 

was  partly  cloudy.  The  two  C^*  indications  are  well  correlated,  and  the 
10.6  micron  curve  responds  to  each  C^*  peak  but  with  a  definite  time  delay 
which  indicates  the  movement  of  general  conditions  toward  the  transmitter. 
Also,  the  4880A  variance,  which  is  saturated,  is  anticorrelated  with  that  for 
10.6  microns,  which  illustrates  the  fall-off  property  "beyond  saturation":  as 
C  2  increases  further,  the  scintillations  actually  decrease.  It  may  be  noted 
that  short-term  turbulence  was  very  large  during  this  run,  with  optical  C^2 
values  of  3  x  10‘1S  m  “2^s  and  10.6  micron  variances  of  as  high  as  0. 15. 


STATISTICAL  STATIONARITY  CONSIDERATIONS  AND  THE 
INERTIAL  SUBRANGE 


There  exist  subtle  relationships  between  statistical  stationarity , 
averaging  times,  low  temporal  or  spatial  frequency  components  of  the 
phenomena,  and  the  interaction  of  the  inertial  subrange  with  the  optical  filter 
function.  Although  the  stationarity  of  atmospheric  phenomena  can  be  treated 

15 

with  some  rigor,  it  is  perhaps  more  fruitful  to  consider  the  direct  physical 
implications.  For  example,  sizeable  fluctuations  of  the  strength  of  turbulence 
(as  averaged  v/ith  e.g.  RC=  100  seconds)  may  result  from  partly  cloudy  skies, 
gusty  winds,  or  less  discernible  causes.  Nonstationarity  may  be  cyclic  or 
monotonic  in  nature.  To  the  extent  that  monotonic  trends  do  not  enter  seriously, 
one  may  use  very  long  averaging  times  to  compare  various  interrelated  experi¬ 
mental  parameters  with  theoretical  predictions.  However,  for  practical 
questions  such  as  the  degree  of  fading  or  the  short-term  error  rates  in  a  com¬ 
munications  or  radar  system,  a  failure  to  follow  or  predict  shorter-term 
extremes  of  turbulence  effects  may  result  in  unrealistic  design  appraisals. 

When  the  medium  is  stationary,  rapidly  fluctuating  phenomena 
such  as  scintillations  at  visible  wavelengths  will  yield  highly  repeatable 
statistical  results  over  e.g.  10  second  averaging  times,  while  phenomena 
such  as  10. 6  micron  scintillations,  with  strong  low  frequency  components, 
require  much  greater  averaging  times  to  reduce  the  experimental  spread. 

It  is  often  somewhat  difficult  to  define  a  difference  between  nonstationarity 
and  large  low-frequency  components,  except  through  the  simultaneous 
monitoring  of  more  rapidly  fluctuating  phenomena. 


*An  experimental  technique  for  obtaining  accurate  probability  sweeps  during 
changeable  conditions  is  to  magnetically  record  the  variable  (AT  or  log  A) 
and  then  measure  discrete  probability  points  by  repeated  scans  of  the  same 
time  window. 
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A  particularly  important  aspect  of  this  question  arises  in 
relating  the  temperature  fluctuations,  which  are  dominated  by  low  frequencies 
and  which  evidence  probability  characteristics  which  might  be  interpreted 
either  as  nonstationarity  or  as  strongly  nongaussian  behavior,  with  optical/IR 
scintillations,  which  are  much  faster  and  usually  log  normal.  It  is  con¬ 
jectured  that  the  nongaussian  probability  tails  for  the  thermal  fluctuations 
are  generally  manifestations  of  (slow)  fluctuations  with  corresponding  spatial 
frequencies  which  are  largely  outside  the  inertial  subrange  and,  in  most  cases, 
outside  the  significant  range  of  the  optical  filter  function.  However,  since 
the  slope  of  the  linear  probability  range  predicts  too  small  a  scintlllation- 

variance,  relative  to  the  short-path  optical  C  2  determination,  it  must  be 

n 

concluded  thai  these  probability  tails  may  not  be  completely  neglected.  A 
direct  determination  of  the  outer  scale  would  be  desirable,  through  large 
thermal  probe  separations  or  very-low-frequency  spectral  analysis  of  the  tem¬ 
perature  fluctuations.  A  related  Item  is  the  optimum  value  of  the  large,  series 
coupling  capacitor  in  the  thermal  probe  instrument,  since  we  are  actually 
measuring  AT  -  <AT>. 

III.  PLANS  FOR  'a  HE  FOLLOWING  PERIODS 

Our  primary  task  now  Is  to  obtain  sufficient  quantities  of  point- 
transmitter  data  under  a  variety  of  conditions  to  substantiate  our  tentative 
conclusions.  Prior  to  the  onset  of  consistent  good  weather,  we  will  attempt 
to  develop  the  thermal  determination  of  the  inner  (and  possibly  outer)  scales 
to  the  point  where  these  parameters  may  become  part  of  a  standard  run. 

After  obtaining  the  large  body  of  fundamental  results,  we  will  obtain  data 
for  varying  transmitter  apertures. 

During  the  late  period  of  the  follow-on  program,  we  may  find  it  feasible 
to  conduct  limited  experiments  over  a  much  longei  path.  These  are  desirable 
in  order  to: 
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(1) 


Demonstrate  10,6  micron  scintillation  behavior  at 
large  variances  (or  saturated) 


and  :(2>  _  .-“Determine  whether  or  not  a  limitation  occurs  in  the 

-•  -  -  value  0i  vs  ,  range  (TO ,•  when  €^0i  becomes  saroxs ted 


(at  ’488oA). 


The  latter  item  is  suggested  by  a  heuristic  argument.  The  fact  that 

r  increases  as  (\  LV*  while  the  "coherence  diameter’1  r  decreases  with 

•  a  .  ,  o 

range,  suggests  that  the  process  of  setting  will  define  a  critical 

'  range  (L^)  fax  any  given  turbulence  level.  As  pointed  out  in  Ref.  4,  the 
'equality  occurs  for  a  particular  C^(0)  «  1  which  is  independent  of  the  values 

of  JL,  C  *,  and  d.  giving  rise  to  this  variance  (This  is  consistent  with  our 

_ •  .  ,n 

new  experimental  evidence  that  the  saturation  level  is  not  waveiength-rieper.dent, ) 
Vft  may  further  speculate  that  r  cannot  continue  to  grow  indefinitely  as  L 
becomes  much  greater  than  L  .  Such  a  limitaDoii  in  r  is  suggested  by  the 
experimentahdata  of  Ref.  T7,  and  has  important  implications  regarding 
aperture  averaging  in  long-path  systems. 

A.  APPENDIX:  NOISE  BANDWIDTH  OF  THERMAL  PROBE 


The  thermal  probe  system  has  been  described  in  detail  in 
recent  reports  on  this  oroject.  Although  the  sensitivity  of  the  unit  is  much 
higher  than  necea sary  for  ground-based  measuremen's ,  the  possible  use  of 
a  similar  system  in  vertical  C. ^  profile  measurements  suggests  that  the 
question  of  ultimate  sensitivity  and  effective  noise  bandwidth  be  discussed 
further . 


The  effective  noise  bandwidth,  or  ability  of  the  instrument  to 
measure  small  relates  to  both  the  response  cpsed  or  i.f.  bandwidth, 

and  to  the  final  integration  time.  If  the  two  bandwidths  are  Bj  and  Bj 


respective!*/,  the  effective  bandwidth  is  basically  propor -iuimi  to  L 


The  question  of  SNR  and  limiting  sensitivity  becomes  a  rather 
subtle  one  because  the  "signal"  is  in  some  sense  the  mean  change  in  the 
meter  reading  (over  system  noise)in  the  presence  of  turbulence,  while  the 
"noise"  relates  to  the  variance  in  this  change  taken  over  an  ensemble.  With¬ 
out  going  into  this  question  in  detail  here,  it  suffices  to  point  out  that  the  use 
of  synchronous  demodulation  is  superior,  in  marginal  detection  situations, 
because  envelope  demodulation  results  in  the  i.f.  noise  suppressing  the  signal 
and  in  unnecessary  baseband  self-convolution  of  the  noise  prior  to  the  squar¬ 
ing  operation. 


In  our  unit,  which  uses  a  low-noise  front-end  technique  but 

envelope  demodulation  for  simplicity,  an  i.f.  SNR  of  unity  occurs  for  a  value 

of  C  2  of  2  x  10~17  m”2/1.  This  is  consistent  with  the  residual  output  (AT)2 
n 

indication  (at  highest  gain)  corresponding  to  system  noise.  If  we  were  to 
approach  this  level  of  C^2  e.g.  at  higher  altitudes,  we  would  want  to  use 
synchronous  demodulation  in  order  to  prevent  the  suppression  of  marginal 
signals  by  the  noise.  In  the  present  case,  this  is  two  orders  of  magnitude 
below  values  which  we  have  seen. 
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FIGURES 


Multiwavelength  scintillations  over  a  one  mile  path.  The 
xor  responding  wavelengths  (top  to  bottom)  are  488GA-1.  lop 
and  48S0A-lG.6u  .  The  horizontal  scale  is  0.1  cm/sec. 

•Simultaneous  probability  distributions  of  log  scintillations  over 
a  one  mile  path  (Run  #4). 

Simultaneous  covariance  results  over  a  one  mile  path  (Run  * 3)  . 

Simultaneous  spectra  of  scintillations  over  a  one  mile  path 
(Run  # 3). 

Examples  of  probability  distributions  of  thermal  fluctuations  (AT). 
Spectrum  of  rms  thermal  fluctuations  (AT)  for  Run  #2. 

Thermal  fluctuations  (Run  ?2).  The  horizontal  scale  is  0.5  sec/cm. 


Log  amplitude  variance  vs.  receiver  aperture  diameter  (4880A). 

Probability  distribution  of  log  amplitude  for  large  receiver  aperture, 
corresponding  to  satured  conditions  on  Figure  8. 

Simultaneous  log  amplitude  variances  (Run  £6). 

Simultaneous  curbulence  strength  measurements  and  log  amplitude 
variances  for  unstable,  large-turbulence  conditions. 
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A  series  of  significant  multiwavelength  laser  propagation  experiments 
have -been.- conducted  with  point-source  transmitters  and  independent  determina¬ 
tions  of  turbulence  strengths.  Under  non-saturat ed  scintillation  conditions, 
M-find  that  the  theoretical  rarelangth  prediction  of  log  amplitude  variance 
is  valid.  Also,  we  find. that  the  near-  and  far-infrared  wavelengths  do  not 
saturate  at  lower  variances  than  .visible,  wavelengths ,  and  large  infrared 
variances  have  been  measured.  Me-  find, that  different  techniques  for 
measuring  the  strength  of  turbulence  yield  different  (out  related)  results. 

..  Jia...have  found  that  the  multiwavelength  covariance  characteristics  are 
usually  independent  of  conditions  and  are  near  theoretical  values,  although 
certain  anomalous  results  suggest  limitations  in  the  extent  of  the  inertial 
subrange  turbulence  model.  -#e^hare  found  that  large  receiver  apertures  do  not 
result  in  the  degree  of  reduction  in  total  signal  fluctuations  that  would  be 
expected  from  the  theory.  1  Einaily,  we  have  found  tnat  simultaneous  records 
of  changing  variances  and  turbulence  strengths  yield  interesting  enaracteristics 
whicn  may  be  related  to  the  6tationarity  of  the  medium. 
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